%% 断裂轨迹补点
function [x,y]=line_point(p1,p2) %p1,p2分别为断裂点，恢复点
x1=p1(1);y1=p1(2);x2=p2(1);y2=p2(2); %断裂点坐标（x1,y1）;恢复点坐标（x2,y2）
a=(y2-y1)/(x2-x1); %参数a
b=p1(2)-a*p1(1); %参数b
if x1<x2
    x=x1:0.0005:x2; 
else
    x=x2:0.0005:x1;
end
y=a*x+b; %构造的直线方程，以及对应的y
end